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1 . Introduction 

Nuclear matter is a hypothetic substance consisting 
of interacting neutrons and protons in large enough 
numbers that the system can be considered infinite. 
The properties of this medium are determined by the 
nucleon-nucleon potential (usually with the Coulomb 
forces between the protons omitted to avoid the in- 
finite Coulomb potential). Two properties, the equi- 
librium density and binding energy per nucleon, can 
be derived from measurements on finite nuclei be- 
cause the conditions at the center of a very large 
nucleus approximate those of nuclear matter. 

The purpose of this investigation is to study the 
important interaction mechanisms of photons in 
nuclear matter as a function of frequency in order to 
estimate the photon propagation characteristics in 
the medium. In the next section we review some of 
the basic facts of nuclear matter. In section 3 the 
photon-nucleon cross sections for various processes 
are computed. These processes are then reexamined 
in section 4 for the case where electrons are added to 
the medium causing net charge neutralization. Finally, 
a discussion of the results and conclusions are pre- 
sented in section 5. 

2. Nuclear Matter 

The charge density distribution of finite nuclei as 
determined by elastic electron scattering indicates 
the nucleon density of heavy nuclei is fairly constant 
in the interior of the nucleus with a value 

p = 0.170 nucleons/fm 3 . 

We shall adopt this value for infinite nuclear matter. 
The total binding energy of finite nuclei as expressed 
in a semiempirical mass formula can be extrapolated 
to an infinite number of nucleons (half protons, half 



neutrons). The average binding energy per nucleon in 
nuclear matter (i.e., the coefficient of the term linear 
in number of nucleons) is 

B = -15.68 MeV/nucleon. 

It is the goal of nuclear matter theory to reproduce 
these "observed" properties, p and B, starting from a 
potential V(r) between two nucleons. A review of these 
calculations is given by Bethe [l] 1 . The Reid soft- 
core potential, one of the better phenomenological 
nucleon-nucleon potentials, yields p RSC = 0.20 fm~ 3 
and # RS( =— 11.25 MeV in a calculation of the two- 
body contributions. Estimates of higher order contri- 
butions seem to account for discrepancies between 
the Reid potential predictions and the observables. 

Nuclear matter is usually discussed within the frame- 
work of the Fermi gas model. For no interactions among 
nucleons (V(r) =0) the individual nucleon wave func- 
tion is a plane wave of momentum k 

These functions are normalized in a unit volume and 
are subject to antisymmetrization for identical nucle- 
ons. Our units are chosen to have H = c=l. 

The distribution of nucleon momentum in the Fermi 
gas model is uniform from k = to A = A> (the Fermi 
momentum). The maximum momentum A> is related 
to the nucleon density by 



k 3 = n 

k F 2 P 

which gives A> =1.36 fm -1 . This value of the Fermi 
momentum is also in agreement with inelastic electron 



1 Figures in brackets indicate the literature references at the end of this paper. 



scattering in the quasielastic region. This value cor- 
responds to a maximum kinetic energy per particle of 

-f- = 37MeV 

2M 

and an average kinetic energy of 22 MeV. These plane 
waves form the basis set used in the solution of the 
problem of interacting nucleons. 

In nuclear matter the Schrodinger equation for the 
wave function i//(r) of two nucleons with relative 
momentum k= ki — 1*2 interacting through a potential 
V(r) is modified by a term which takes into account 
the fact that nucleons cannot scatter into momentum 
states already occupied by other nucleons. The modi- 
fication as derived by Bethe and Goldstone for a 
s-wave potential is 



dr 2 



+ k 2 u(r)=MV(r) u(r) 
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where u(r) = kr t//(r) and 

sin /^(r — r') sin &/, (r + r') 



X(r,r') = 



r+r' 



The solution of the Bethe-Goldstone equation can be 
written as a free wave minus a defect part 

u(r) = sin kr-ku(r). 

Figure 2 shows the defect function [2] for the Reid 
soft-core potential. The function Aw/ k is not sensitive 
to either the relative momentum, k, or the total momen- 
tum ki + k 2 . An important feature of Au(r) is that it is 
small by the time r approaches the average inter- 
nucleon separation distance, p~ 1/3 = 1.8 fm. This 
means that for subsequent collisions nucleons approach 
each other in relatively pure plane-wave states. 
The defect function normalization is such that 

K = -pj X \Au(r)\*dr 

is the probability that a s-wave pair of particles is 
outside the Fermi sea. For the Reid soft-core potential 
k= 0.14 summed over all partial waves. 

The Pauli exclusion principle which prevented 
nucleon-nucleon scattering to occupied momentum 
states will also modify photon-nucleon interactions. 
If in free space a photon process imparts a momentum 
transfer q to a nucleon, the cross section in nuclear 
matter will be reduced by a factor 



R(q) = 



4kr 



1 
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because only a fraction R of the nucleons will have an 
initial momentum ksuch that the final momentum will 
exceed the Fermi momentum | k + q | > k F . The 
expression is derived by calculating the excluded 
volume of the initial and final Fermi momentum 
spheres shown in figure 1. 




FIGURE 1. Initial (left sphere) and final (right sphere) momentum 
distributions in a Fermi gas model in which a momentum q is to 
be given each nucleon. 

Only those nucleons for which | k+ q | > ky are allowed to make the transition. The 
fraction R(q) is the ratio of excluded volume (shaded) to the total (4/37T k F 3 ). 

3. Photon Interactions 

Consider a photon of frequency co incident on a 
semi-infinite slab of nuclear matter. The index of 
refraction (ratio of wavelength in vacuum to that in 
the material) is related to the forward elastic scattering 
amplitude /( co ) by 



n 2 (to) = l + 



47rp/(w) 



q < 2k F 



In general, both / and n are complex numbers. The 
imaginary part of / is related to the absorption cross 
section through the optical theorem 

47T Im/(co) = cocr(co). 

The real part of /can be calculated from the imaginary 
part by the use of a dispersion relation, but we will 
not make use of this approach in our estimates. We 
now examine the photon interaction mechanisms [3]. 

3.1. Thomson Scattering 

Electromagnetic waves scatter from the free proton 
(mass M, charge e) with an amplitude 



/=-- = -L54X10- 3 fm, 
M 
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which is independent of frequency. The forward ampli- 
tude is the same in nuclear matter because no momen- 
tum is transferred to the proton. It is convenient to 
introduce a quantity called the plasma frequency 

e 2 

Col = 4<7Tpp — > 

which for p P = 0.085 protons/fm 3 has the value coo = 
8.0 MeV. The index of refraction due to Thomson 
scattering can be written as 

col 

!*(«) = 1 "°. 

co- 
Note that for co < 8 MeV the index n is imaginary 
which means an incident photon will be reflected as it 
tries to enter the medium. This phenomenon is analo- 
gous to the reflection of light waves from the surface 
of a metal due to Thomson scattering by the conduction 
electrons. 

Within the medium the photon intensity falls off 
exponentially and the 1/e penetration depth is given by 



d T \ 



i- (^o 



-1/2 



For co < < coq this distance is 25 fm. 

Above the plasma frequency n{co) is real and the 
photons propagate freely with the index approaching 
unity from below. 

3.2. Proton Compton Effect 

Photon scattering by the proton's charge is inelastic 
at angles other than 0° and therefore absorptive. The 
differential cross section for scattering by a free 
proton is given by the Klein-Nishina formula 



J free / p 2 \ 2 ] 
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where the scattered photon energy is 
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For incident photon energies small compared to the 
proton mass, co < < M =938 MeV, the photon energy 
loss co — co' is small and the free cross section can be 
approximated by 

J free / p i x 2 1 

uu Compton f \ /i . 2D\ 

— 77: — —It; -(1 + cos 2 0). 
dn \ M I 2 

In nuclear matter this cross section must be multi- 
plied by the Pauli reduction factor. The momentum 
transfer to the proton is 

1 
q = 2co sin - 6. 



The total Compton cross section for a noninteracting 
Fermi gas is then 
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For example, cr™ t (100 MeV) = 0.73 X lO" 31 cm 2 . 

r Compton / 

3.3. Pair Production 

An important absorption mechanism is the conver- 
sion of the photon to an electron-positron pair in the 
Coulomb field of the proton. The free proton cross 
section for co > > m e is given by 






28 2co 218 
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where m e is the mass of the electron. The proton is 
required to take up an average momentum transfer of 
q = m e . Thus, the nuclear matter cross section per 
proton is 



^NM r 1 _^ Z?Ie free, . 

(J (co) — CJ (co) , 

pair 4 h F pair 



A 1/e absorption length can be computed using 

d= (prcr)- 1 . 
Typical values are: 



CO 

(MeV) 


NM 

u pair 

(Xl0- 27 cm 2 ) 


Imfpair 

(Xl0- 3 fm) 


"pair 

(fm) 
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30. 

100. 


0.38 

0.77 
1.19 


0.15 
0.93 
4.80 


306. 

153. 

99. 



The imaginary part of the forward scattering ampli- 
tude exceeds the real Thomson aplitude above 50 
MeV. 



3.4. Absorption on Correlated Neutron-Proton Pairs 

Single particle ejection, which is an important 
photoabsorption mechanism in finite nuclei, is not 
possible in a noninteracting system in which the 
nucleon has definite linear momentum because energy 
and momentum cannot be simultaneously conserved. 
However, in an interacting Fermi gas the defect 
wave function has a spread of momentum components 
that allows two nucleons to absorb a photon and be 
ejected from the Fermi sea. 

An estimate of the absorption cross section of this 
process can be made by considering the electric 
dipole transition of a correlated neutron-proton 
pair using the s-wave defect wave function as the 
initial state and a free p-wave for the final state 
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The density of final two-nucleon states is a function 
of the final relative momentum p 



cPp 



dn 



do) (27r) 3 dco (2tt) 3 2 
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The El matrix element is calculated from the momen- 
tum space wave functions as 



-/ 



where e is the photon polarization vector. 
A properly normalized initial state is 

fc(A) = (3p/8V*-A*(A) 

where for convenience we use the Fourier transform 
of the defect wave function 



ty(k) 



(2tt) 3 / 2 J 



• Au(r) 

" k r — d 3 k. 
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The final state with both particles outside the Fermi 
momentum sphere is 

l\f f = (27r) 3 / 2 8 3 (p-k). 

The matrix element is simply evaluated as 

M{p) = (3 p/8) 1 / 2 (2tt) 3 / 2 € • pAi//(p). 

The total cross section averaged over photon 
polarization is then 



l \- 3^ (P 3 



W(P) 



We take the threshold for the reaction to be twice the 
average single particle binding energy, 16 MeV. 
Thus, photon energy and final relative momentum 
are related by 



a> = ^ + 32 MeV. 



The resulting cross section is plotted in figure 2. At 
the peak co = 100 MeV the nuclear matter cross section 
is about 13 percent of the free deuteron cross section. 

3.5. Pion Production 

At photon energies greater than the mass of the pion 
o) > m n = 140 MeV the reaction 
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FIGURE 2. Upper right insert shows s-wave defect wave function for 
a soft-core nucleon-nucleon potential. 

The curve cr pn is the cross section per nucleon for the absorption of a photon by a 
correlated neutron-proton pair calculated using Au. 

occurs for both neutrons and protons. The cross section 
for free nucleons reaches a peak at co = 300 MeV of 
0.5 X 10 -27 cm 2 (the 3, 3 resonance). After a few 
smaller resonances the cross section is essentially 
constant between 2 GeV and 20 GeV at 0.1 X 10 " 27 
cm 2 . One might expect that the nuclear matter cross 
section would be the free cross section modified by 
the Pauli reduction factor R(q) where the nucleon 
recoil momentum is 



Im^M 



M + 



(c 



-11/2 



However, the A dependence of the high energy ab- 
sorption cross section for finite nuclei 



cry a 
cryN 



A - 75 , co = 0.14 to 2 GeV 
A™\ co > 2 GeV 



y + N -> N+ir 



indicates [4] a shadowing phenomenon which is more 
complicated than the incoherent sum of individual 
nucleon events in nuclear matter. Present theories 
contend that a high energy photon spends a fraction 
of the time as a spin-one meson (vector dominance 
model). Mesons propagating in nuclear matter have 
a large and complex index of refraction [5]. For our 
purposes, however, it is probably safe to assume that 
the photo pion absorption cross section per nucleon 
is no larger than that of the free nucleon at the same 
photon energy. In this case, the high energy absorption 
will be dominated by the ever-increasing e + e~ pro- 
duction process. 

4. Photon Interactions in Charge Neutralized Nuclear 
Matter 

If electrons are added to the proton-neutron matter 
to achieve charge neutrality (p e = pp — 0.085 fm~ 3 ) 
the optical properties are drastically changed. Miiller 
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and Rafelski [6] argue that charge neutralization occurs 
spontaneously when the Coulomb potential exceeds 
[k F 2 + m e yi 2 - k F = 268 MeV. This occurs for finite 
nuclei which have more than 10 4 protons. 

The plasma frequency of the electrons is much 
higher than that of the protons because of the smaller 
electron mass, 

^electron = (^ p<? ^j^ ) 1/2 = 342 MeV, 

with an associated penetration depth less than a fm. 
As a consequence, any photons with energy less than 
342 MeV will be reflected from the surface of a semi- 
infinite slab of this material. 

Photon propagation in the medium becomes possible 
again for a> > 342 MeV. Pair production in the Coulomb 
fields of the electrons and protons is influenced by 
the requirement that the e~ which is produced must 
have a momentum greater than A>. Thus, the threshold 
for pair production is raised from 2m e to Ay and the 

In — term in the high energy cross section formula 
m e 

becomes In - — . In the threshold region the e + and e~ 

energies are not shared symmetrically and there will 
be mutual screening of the electron and proton 
Coulomb fields. An estimate of the combined pair 
production cross section can only be made for co > 2 
GeV, viz., 
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The photo pion absorption cross section is unaf- 
fected by the electron component of nuclear matter. 
Above 2 GeV the pair cross section is larger than 
C777T, but the reverse is probably true between 342 
MeV and 2 GeV. Compton scattering on the electrons 
is smaller than both these cross sections. 



nuclei occurs coherently in the intense Coulomb 
field at the edge of the nucleus rather than from 
individual protons in the interior. 

Electron-proton-neutron matter has a much higher 
plasma frequency so that photons do not penetrate 
the medium until a> > 342 MeV. Pion production from 
the protons and neutrons dominates the absorption 
up to about 2 GeV, then pair production becomes more 
important; however, both these processes are difficult 
to calculate accurately. 
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FIGURE 3. Cross sections for dominant photoabsorption processes 
in nuclear matter. 

The vertical line at w = 8 MeV marks the frequency below which photons will not prop- 
agate in matter composed of protons and neutrons. The line at o> =342 MeV is the limiting 
energy for matter composed of electrons, protons, and neutrons. The solid curve laheled 
a- p a i r is the absorption cross section per proton for e e~ production in pn matter: the dashed 
curve for epn matter. The curve labeled (Ty n is the free nucleon photo pion absorption cross 
section. This cross section is probably an upper limit to the process in nuclear matter 
Other photon interaction cross sections, Compton scattering, and nucleon pair absorption, 
are less that 10~ 2! ' cm 2 at all energies. 
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5. Discussion and Conclusions 

The main results are displayed in figure 3. Proton- 
neutron matter is dominated by reflective Thomson 
scattering below 8 MeV and pair production above the 
proton plasma frequency. All real nuclei are too small 
to exhibit the low frequency imaginary index of re- 
fraction expected of np matter since the largest 
diameter (14 fm) is still small compared to the pene- 
tration depth (25 fm). Also, pair production on finite 
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